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ABSTRACT 



A star can be tidally disrupted around a massive black hole. It has been known 

CS ' that the debris forms a precessing stream, which may collide with itself. The stream 

collision is a key process determining the subsequent evolution of the stellar debris: if the 

Q\ . orbital energy is efficiently dissipated, the debris will eventually form a circular disk (or 

torus). In this paper, we have numerically studied such stream collision resulting from 

Ph' the encounter between a 10^ Mq black hole and a 1 Mq normal star with a pericenter 

Q ■ radius of 100 Rq. A simple treatment for radiative cooling has been adopted for both 

optically thick and thin regions. We have found that approximately 10 to 15% of the 

initial kinetic energy of the streams is converted into thermal energy during the collision. 

The spread in angular momentum of the incoming stream is increased by a factor of 2 to 

3, and such increase, together with the decrease in kinetic energy, significantly helps the 

circularization process. Initial luminosity burst due to the collision may reach as high as 

_c^ \ 10^^ erg sec~^ in 10*^ sec, after which the luminosity increases again (but slowly this time) 

to a steady value of a few lO^^ergsec"^ in a few times of 10^ sec. The radiation from 

the system is expected to be close to Planckian with effective temperature of ~ 10^ K. 

Subject headings: black hole physics — hydrodynamics — radiation mechanisms: ther- 
mal — methods: numerical — galaxies: nuclei 
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1. INTRODUCTION 

There are dynamical evidences for supermassive 
black holes in nearby galaxies including our own (e.g., 
Eckart & Genzel 1996). Since these black holes are 
embedded in a dense star cluster, tidal disruption of 
stars can occur at a rate determined by the stellar 
density and velocity dispersion. Rough estimates in- 
dicate around 10~^ per year for disruptive encounters 
between a star and a massive black hole for the Milky 
way (e.g., Goodman & Lee 1989). 

The debris from the disruptive encounter forms a 
precessing stream, and self-interaction (collision be- 
tween different parts of the stream) may occur soon 
after the most tightly bound debris begins its second 
orbit. Evans & Kochanek (1989), Kochanek (1993), 
Laguna et al. (1993), and Monaghan & Lee (1994) 
have tracked the debris from the disruption to the 
collision phase. The disruption can in principle lead 
to the flare up of luminosity in galactic nuclei (e.g., 
Rees 1988) followed by an accretion phase lasting 
for several years (Cannizzo, Lee & Goodman 1990; 
Loeb & Ulmer 1997; Ulmer, Paczyhski, & Goodman 
1998). The studies on this rather long-lasting, ac- 
creting phase started from the time when the debris 
is circularized assuming that such circularization oc- 
curs in a fairly short amount of time (several times 
of the orbital period of the most tightly bound de- 
bris). The present study targets the evolution of the 
debris during and right after the collision. The circu- 
larization process will take much longer than the time 
scale that we study here, and is beyond the scope of 
the present study. 

The detailed evolution of the stellar debris has not 
been studied carefully because of the intrinsic com- 
plexity of the problem. The difhculty stems from the 
huge dynamic range of the physical parameters of the 
debris around the black hole. For example, the peri- 
central distance of the debris' orbit is typically 1 AU 
while the apocenter distance could be of order of pc. 
The most important phenomenon in determining the 
fate of stellar debris is the collision between the in- 
coming and outgoing streams due to relativistic pre- 
cession (e.g. Rees 1988, Cannizzo et al. 1990). The 
supersonic collision between two streams can in prin- 
ciple convert orbital energy into thermal energy. If 
this process efficiently removes the orbital energy, the 
stream transforms itself into a circular disk (or torus) 
around the black hole. The stream crossing is likely 
to be intermittent in nature because the collision can 



destroy the stream geometry (Lee, Kang, & Ryu 1995; 
LKR hereafter) . Such a cycle will repeat with period 
of free-fall time at the crossing point. 

LKR have numerically studied the stream collisions 
and showed that the shock should be able to convert 
the ordered orbital motion of two fluids into expansion 
of a bubble and thus produce the torus around the 
black hole. For simplicity, LKR assumed that the col- 
lision is adiabatic and gas pressure is dominant over 
radiation pressure. However, it is more likely that the 
radiation pressure is more important at the shocked 
region and radiative cooling may play an important 
role. The purpose of the present study is to extend 
LKR's numerical simulations including radiative pro- 
cesses. We pay special attention to the physical status 
of shocked region and expanding post-collision gas, 
the amounts of thermalization and radiative cooling, 
and changes in orbital energy and angular momentum 
of the streams. 

We describe our numerical methods for hydrody- 
namics and radiative cooling in §y. Simulation results 
and energy conversion issue are discussed in ^ and 
some estimates for orbital and radiative evolution are 
made in ^. The final section summarizes our find- 
ings. 

2. NUMERICAL METHODS 
2.1. Hydrodynamic Simulation 

We perform simulations of supersonic collisions be- 
tween two gas streams with a three-dimensional hy- 
drodynamics code based on the total variation di- 
minishing (TVD) scheme (Harten 1983; Ryu et al. 
1993), which is an explicit, second-order, Eulerian 
finite-differencing scheme. The code used in LKR has 
been modified here to include the effects of radiative 
cooling. The treatment for the radiative cooling is 
discussed in the next subsection. 

While hydrodynamic simulations in LKR were per- 
formed purely dimensionlessly, the introduction of 
radiative processes requires some dimensional con- 
straints. In this paper, all simulation variables have 
specific physical values. We targeted the stream de- 
bris resulting from the encounter between a 10^ Mq 
black hole and a normal star with mass of 1 Mq and 
radius of 1 Rq, which initially has a parabolic orbit 
with a pericenter radius of 100 Rq. 

Some initial conditions of our simulations were ob- 
tained from the results by Kochanek (1994) and Lee 



& Kim (1996). Velocities and mass fluxes of two 
streams were obtained from the values at i « 1750 to 
of his Figure 6 {vi = 0.0147 c, V2 = 0.0083 c, Mi = 
lM0yr~^, and M2 — O.5A/0yr^^). We have chosen 
this epoch because AI2 becomes significant compared 
to Ml (M2 ~ Ml/2) at this epoch (subscript 1 is 
for the stream making its first orbit, stream 1, and 2 
is for the stream making its second orbit, stream 2). 
Kochanek found that the angle between the streams 
at the collision point is typically 130° to 140° and 
here we adopted a collision angle of 135°. However, 
his results for the widths of two streams Ai and A2 
at a certain epoch depend on the way of treating the 
viscosity. For this reason, we perform simulations for 
five different sets of A's. We have fixed the heights of 
two streams. Hi and H2 except for one run because 
these quantities are less dependent on the viscosity 
treatment. 

Model parameters of our simulations are shown in 
Table |[ Our standard model is run 7, where the 
cross sections of the streams are round (Ai = Hi & 
A2 = H2), two streams move in the same z = 50 plane 
(hereafter the coordinate positions are represented by 
the pixel number; see Iceii below for the physical size 
of each pixel), and the mass flux and the radius of 
stream 1 are 2 and 5/3 times larger than stream 2, 
respectively. Runs 2, 8, 9, and 10 are for different 
A and H values. Runs 11 and 12, where stream 2 
is set to move in z = 54 and 58, respectively, are for 
the estimation of the effects of off-centered collision 
that could occur when there is Lense-Thirring preces- 
sion. Finally, to see the effects of radiation, we will 
compare runs 7 and 1, where radiative cooling is not 
considered. 

The computational box of our simulations consists 
of 101'^ cubical cells. The physical length of each cell 
Iceii is 5.93 X 10^° cm. The gas is assumed to be ionized 
already before the collision and the initial pressure of 
the streams is set such that the Mach number of the 
stream is 300. This gives an initial temperature of 
about 2 X 10'* K for stream 1 and 7 x 10'^ K for stream 
2. The initial density of the ambient medium is 10~* 
times smaller than the average density of the streams 
and the initial temperature of the ambient medium 
is set to be slightly smaller than that of the streams. 
The initial density profile of stream's cross section has 
a functional form of exp(— 4y^/A^ — 4z^/H'^). Unlike 
the simulations of LKR, where two streams are set 
to collide on their front ends, our stream 2 is set to 
collide on the side of stream 1 , which more resembles 



the geometry of actual stream collisions. Streams 1 
and 2 are set to keep flowing into the simulation box 
until the end of the simulation. 

The total integration time of our simulations is 
4 X lO'* sec. Kochanek (1994) finds that for the en- 
counter with parameters adopted in this paper, the 
debris stream will self-intersect at 5000-6000 i?© from 
the black hole on its second orbit and the first col- 
lision phase will last for ^ 10^ sec. Thus the one- 
dimensional size of our simulation box is about 100 
times smaller than the distance to the black hole from 
the collision region and our simulations cover only 
small fraction of the first collision phase. 

2.2. Radiative Processes 

The average particle number density of the streams 
is of order 10^^ cm"'^ and the characteristic radius of 
the streams is of order 10"'^^ cm. For the typical col- 
umn density of ^ lO^^cm"^, the electron scattering 
optical depth Tes is much larger than unity, and so is 
the effective optical depths r* = [3t//(t// -|- Tes)]^/^ 
(Rybicki & Lightman, 1979), where t// is the Rosse- 
land mean free-free absortion optical depth. Thus the 
post-collision region may be approximated by thermal 
equilibrium where matter and radiation arc strongly 
coupled. The total internal energy density, u, is now 
the sum of the gas thermal energy density, Ug, and 
the radiation energy density, Ur'. 
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where n is the gas number density, k the Boltzmann 
constant, a the radiation constant, and 7g (=5/3) the 
adiabatic index for gas. The mixture of gas and radi- 
ation gives an effective adiabatic index 7e// between 
5/3 and 4/3, depending on the ratio Ug/ur (or on the 
ratio Pg/Pr, where Pg and Pr are the pressure of the 
gas and the radiation respectively). 

We assume that the post-collision region has phys- 
ical conditions P^ ^ Pg and adopt j^ff — 4/3 ex- 
cept for run 1. The validity of this assumption may 
be checked in Figure |l], which illustrates the density 
weighted distribution of post-collision 7e// for run 7 
at i := 4 X 10* sec. Here, 7e// is defined by 
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where Pg = nkT and Pr = ^aT"*. Most of the ther- 
malized material falls onto the region where "feff is 
close to 4/3. 



Since the generic temperature of the post-collision 
region is well above 10'* K, we assume that gas is 
completely ionized and is cooled by bremsstrahlung 
when r* < 1. Between 10^ K and 10^ K, line cooling 
by various atoms could be significant compared to 
bremsstrahlung. However, the opacity for line pho- 
tons is usually much larger than that for the contin- 
uum photons and therefore do not contribute much 
to actual cooling. 

Bremsstrahlung photons can not escape freely ei- 
ther in most cases. They have to diffuse out from the 
inner, generally hotter, region to the outer, generally 
cooler, region, and in doing so heat up the outer parts. 
As in the stellar interior, the amount of diffusion de- 
pends on the temperature gradient, and the temper- 
ature gradient itself is maintained by the amount of 
energy transported from one part of the gas to the 
other. However, one important difference between 
colliding gas streams and the stellar interior is that 
part of the bulk kinetic energy can always be con- 
verted to the thermal energy in any part of the flow 
which is ultimately radiated away, whereas only the 
nuclear burning core generates the energy and outer 
envelope passively transfers the radiation in the stel- 
lar interior. Hence, the correct way to handle this 
kind of dynamic mixture of gas plus radiation is to 
fully solve the three-dimensional radiative transfer at 
a given instant along with hydrodynamics. This will 
tell us how much radiation is generated and trans- 
ported from one part to the other and, therefore, de- 
termine the physical state of the gas and radiation in 
the next time step. However, this is quite a formidable 
task and we need to find a simpler way to treat this. 

One reasonable and efficient way for this type of 
hydrodynamic simulation is to use a volume cooling 
rate which can approximate the radiation transport 
process. We adopt the following form for the cooling 
(Liang & Wandel 1991; Wandel & Liang 1991): 



e = 7?/3 



rr:4 



1 - e-^' 



R, e-^* +(l~e-^-)77r,/(12CT)' 



(3) 



where 



^-^ = 8.2x10-^; p. 



"// 



i{aff+aes)' 



(4) 



Here e is the cooling rate per unit volume, a the 
Stefan-Boltzmann constant, qb the frequency average 
of the velocity averaged free-free Gaunt factor, and 
gn the Rosseland mean free-free Gaunt factor. Also, 



a// and aes are the Rosseland mean free-free absorp- 
tion and electron scattering coefficients, respectively. 
This form is valid both for effectively optically thin 
(t* <C 1) and for effectively optically thick (r» ^ 1) 
cases. 

For effectively optically thin cases, equation p) be- 
comes the usual brensstrahlung emission rate per unit 
volume, 



e// 



for T* ^ 1- 



(5) 



For effectively optically thick cases, equation (H) re- 
duces to the diffusive cooling rate by thcrmalizcd gas. 
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for r* > 1, 
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where the total optical depth Tfot = Tff + t^^s and 
Rtot = Ttot/{oLff + aes)- Integration of e over the 
whole volume yields roughly the surface area times 
(jT^ff, where T^tr — T'^Ttot- This is equivalent to 
finding the effective temperature of the photosphere 
using the Eddington approximation, and summing up 
the flux over the whole surface. 

The amount of cooling is calculated explicitly by 
equation (^), and is subtracted from the total inter- 
nal energy of each cell after every hydrodynamic step. 
To obtain optical depths, appropriate absorption or 
scattering coefficients were integrated along ±a;, ±2; 
directions, and the smallest among the integrated val- 
ues is chosen. The effective photon travel length R^, 
is determined by T,/[3a//(a// -|-aes)]*^^, where local 
values were used for the calculation of a/y and a^s- 

The exact time-dependence of escaping photons 
from optically thick medium was solved by Sunyaev & 
Titarchuk (1980) for a spherical plasma cloud. Their 
work shows that when a point source, located at the 
center of the cloud of radius R and the optical depth 
T, emits flash of radiation at time t = 0, the photon 
escape rate is highest when t ~ t^sc and exponen- 
tially declines afterward with the same timescale igscj 
where tf-sc — (3/7r^)trf and td = {R/c)t is the pho- 
ton diffusion timescale. However, when the sources 
of photons are homogeneously distributed within the 
sphere, photon escape rate is highest at i = and de- 
clines afterward with the same exponential behaviour. 
When the distribution of photon sources follows a sine 
curve, the exact solution of the escape probability is 
found to be exp(— i/iesc)- Therefore, the fraction of 
photons that escaped from i = to At (^ t^sc) is 
(7r^/3)[At/(Ti?/c)]. This is equivalent to equation 



(||) with a different constant factor. Since sources of 
photons in our simulation are distributed but do not 
foUow the exact sine form, we wiU set this constant 
factor as unity in our simulations. 

Energy of emitted photons can be changed via 
Compton scattering since electron scattering optical 
depth is usually quite high. However, this happens 
only when the typical length scale for absorption is 
longer than that for Compton upscatter. Otherwise, 
photons will be absorbed and thermalized before be- 
ing significantly upscattered. Since the mean free 
path due to electron scattering is (71^(7^3)'^ and pho- 
tons need to be scattered several Am^c^ jkT times to 
gain significant energy boost, the typical length scale 
for upscattering is Acompt = ("■eO'es4fcT/meC^)~^. 
Here, rig is the electron number density, tJos the 
Thomson cross section, and rrLf, the mass of electron. 
In the thermalized regions of our simulations, the ab- 
sorption length scale Aabs — ctjr is much larger than 
Acompt and Compton scattering can be neglected. 

It is possible in high temperature plasma that ions 
decouple from electrons and each develops widely dif- 
ferent temperature. If so, the gas would have to 
be treated as such. During the collision, ions and 
electrons may reach the highest possible tempera- 
ture, which is close to the virial temperature: for 
ions, T™°^ = rripv'^/Sk, and for electrons, T""»^ = 
meV^/3k, where rup and me are the proton and elec- 
tron mass, respectively, and v the collision velocity. 
With our initial stream parameters, we find that the 
Coulomb time scale between ions and electrons at 
virial temperature, ~ 10~^ sec (Park 1990), is a few 
orders of magnitude smaller than the bremsstrahlung 
cooling time scale, ~ 10~^ sec. Besides, the collision 
region is effectively optically thick at this electron 
temperature, and ions will be stronly coupled with 
electrons and radiation will be locked with electrons 
at all times. 



SIMULATION RESULTS 



3.1 



Morphology and Structure 

Figure || shows n (contour) , T (greyscale) , and Vxy 
(velocity projected on the x-y plane; arrow) in z=50 
plane of run 7 at t — 4 x 10^ sec. Stream 1 coming 
from the top boundary is colliding near x,y = (55, 60) 
with stream 2 coming from the lower-right corner. 
Mass elements are almost instantly thermalized in 
the collision area and are spread out in two direc- 
tions. While part of stream 2 is bounced by stream 1 



and heads for -l-a;, -\-y direction, the rest merges with 
stream 1, which is deflected toward —x direction after 
the collision. 

The thermalized gas emerges from the narrow, 
slab-like shock region, and expands out into a larger 
volume in two opposite directions due to the post- 
collision pressure. Thus the post-collision material 
forms two expanding streams several times thicker 
than the pre-collision streams. Hereafter, the post- 
collision gas (PCG) in the -y direction from the colli- 
sion region will be called PCG A, and the other one 
in the opposite direction will be called PCG B. 

Profiles of n, T (defined as eq. O), radiative cool- 
ing per cell (A), Ttot, and t, in z=50 plane of run 7 
at i = 4 X 10'' sec are plotted in Figure for four 
different values of y. Since we did not consider the 
cooling under T = 10"* K, T is not allowed to become 
below lO'' K. Also, T scarecely goes higher than 10^ K 
because cooling becomes more efficient at higher tem- 
perature. On the other hand, r* is greater than 10 in 
most regions because of the high Thomson opacity, 
and goes up as high as ^ 10^ in the regions that are 
not thermalized yet. 

Figure ||(c) shows the profile of y—55, z—50 line 
which includes the central region of the collision. The 
shock front is formed at 40 < x < 50, where temper- 
ature is one to two orders of magnitude higher than 
that of pre-collision gas. Cooling is most prominent 
at the shock slab, but the ratio of the cooling to the 
total internal energy (gas -I- radiation) in that region 
during one simulation time step is only ^ 0.01 %. 
Therefore the radiative cooling time scale is consider- 
ably longer than the hydrodynamical time scale (one 
simulation time step). 

Kochanek (1994) found that the temperature of 
the streams right before the collision is well below 
10"* K. Since we do not consider cooling processes un- 
der 10* K in our simulations, no cooling is expected 
in the pre-collision gas. However, pre-collision stream 
regions, such as 55 < x < 80 of (b) and a; > 80 of 
panel (a) of Figure 0, show some cooling although 
very small. This is because the initial temperature 
of the streams is ~ 2 x lO^'K, which is determined 
by setting Mach number equal to 300. The stream 
material of this temperature will already begin to ex- 
perience radiative cooling even before the collision, 
which is not desirable in our simulations. Higher 
Mach number would give lower initial stream tem- 
perature, but numerical difficulties prohibit us to use 
such high Mach numbers. In any case, the amount 



of such undesirable cooling is negligible compared to 
the one in the pre-collision gas. 

At the boundaries between the expanding PCGs 
and cold streams, such as at a; « 80 of panel (a), 
X ~ 85 of panel (b), and x w 30 of panel (d), the 
flows of the PCG and the stream move in opposite 
directions and shear viscosity thermalizes the bound- 
ary region between the flows. Since such regions have 
relatively small optical depths, thermalized energy is 
efficiently radiated away and local cooling maxima 
occur in those regions. 

A grey scale map of cooling in z=50 plane for the 
same simulation at the same epoch as Figure || is plot- 
ted in Figure H along with the density contour map. 
Cooling is most efficient at the collision region and at 
the mouths of two PCGs. Although our simulations 
do not allow the exchange of photons between gas el- 
ements, we may anticipate that the photons created 
in the collision region will preferentially escape into 
the PCGs rather than into the incoming cold streams 
because the expanding PCGs have much lower den- 
sity and thus play the role of window to the escaping 
photons. The leakage of the thermal energy in the col- 
lision region into the cold, incoming streams will be 
inefficient because the cold streams have huge optical 
depths. 

Profiles along the line connecting points Pi and 
P2 in Z—5Q plane of Figure ^ are shown in Figure |5|, 
where r is the distance in Iceii from point Pi. Note 
that the density profile shows an exponential behav- 
ior at r > 20. Since T and A are tightly related 
to n, these variables also decrease exponentially in 
the same region. The profiles of these variables are 
found to be well fitted by n oc exp(— r/1.7 x 10^^ cm), 
T oc exp(— r/4.9 x 10^^ cm), and A oc exp(— r/9.2 x 
10" cm). 

As we will discuss in the following subsection, 
the fraction of the energy radiated in the simulation 
box to the initial kinetic energy is only ~ 0.3 %. 
Thus the dependence of T on n will not be far from 



the adiabatic relation, T (x n'^ ^ 



(X n 



1/3 



(recall 



that we adopted 7 = 4/3 for all our simulations). 
Indeed, the T scale length is almost one third of 
the n scale length. If the PCG is assumed to ex- 
pand in two dimensions (as when a stream with in- 
finite length expands to the width and height direc- 
tions), one would have Rtot oc n~^/^. Then the A 
scale length, 9.2 x 10^"'^ cm, is a consequence of the 
cooling formula for optically thick region, equation 



(i): A oc T^Rtol^ol^ r4i?-2a-//V T^Rt^ln-'^' 
ex exp(— r/9.0 x 10"'^^ cm), where t^s ^ Tff has been 
assumed. 

Simulations representing two opposite A1/A2 cases, 
runs 8 and 9, show very distinct morphological evo- 
lution as shown in Figures |6| and 0. First, while the 
structure of the thermalized slab in the collision re- 
gion of run 8 is considerably bent toward stream 2, 
that of run 9 is nearly straight with only a small cur- 
vature toward stream 1. Such a difference is a con- 
sequence of different mass flux per unit area, M / dA. 
In run 8, stream 2 has a greater M /dA than that 
of stream 1, thus it penetrates into stream 1 more 
deeply than in run 9 and forms a shock front almost 
perpendicular to both streams although the greater 
M of stream 1 determines the eventual direction of 
PCG. In run 9, on the other hand, M /dA of stream 
2 is significantly smaller than that of stream 1 , there- 
fore stream 2 is effectively bounced at the long shock 
front. The expansion aspects of their PCGs are ac- 
cordingly different. The PCG A of run 8 is largely 
extended, but that of run 9 maintains its slab-like 
shape until it escapes from the simulation box. Thus 
the orbital evolution of the PCG A is expected to be 
considerably different for two simulations, and will be 
discussed again in §^. 

If the central black hole is rotating, the orbital 
plane of the stream could experience the Lens-Thirring 
precession that causes slightly off-centered collision. 
Figure || shows an isodensity plot at i = 4 x lO'^ sec 
for run 12, where only the lower part of stream 2 is 
set to collide with the upper part of stream 1. The 
collision produces an extended PCG A and a narrow 
PCG B with a pitch angle of < 45° from x-y plane. 

3.2. Energy Conversion 

A certain fraction of stream's kinetic energy is con- 
verted into thermal energy due to supersonic collision 
and the key interest of the present study is how much 
fraction of such thermal energy escapes through ra- 
diative cooling. If the collision significantly alters the 
orbital motion of stream 1, the collision phenomena 
becomes intermittent, and the first collision would 
last for about 10^ sec. However, since the hydrody- 
namic simulations in this study are limited within a 
small region around the collision point, the duration 
of simulations with reasonable accuracy is only order 
of the time scale for the stream to cross the computa- 
tional domain, which is few percents of the duration 



of the first collision phase. For this reason, we per- 
form our simulations until energy conversion enters 
a quasi steady state, and extrapolate the results to 
make some estimates on the energy budget problem 
during one whole collision phase. 

Since mass elements are continuously coming into 
and going out of the calculation volume, and since 
it is very difficult to keep track of mass elements in 
the simulations based on fixed data points (grids), 
calculating the amount of energy conversion of a 
mass element which experiences collision is not triv- 
ial. Therefore, we estimate the energy budget in a 
time-cumulative way: Ein^sum is defined as the cu- 
mulated (from t — 0) input energy {Ein^gy^^n = at 
t — 0), -Efh,i«m is the summation of the total internal 
energy (gas-|-radiation) inside the calculation box at 
a certain epoch and cumulated total internal energy 
which has escaped from the box until that epoch, and 
Er,sum is the cumulative energy which has been ra- 
diated away inside the simulation box until the same 
epoch. 

The total internal energy of a mass element before 
collision is negligible compared to its kinetic energy. 
Thus the ratio Rth = {Eth,sum + Er^sum) / Ein,sum is 
the fraction of total thermal energy converted from 
the total input kinetic energy. And the ratio Rr = 
Er,sum/ Ein^sum is the fraction of total radiated en- 
ergy inside the simulation box out of the total input 
energy. The evolution of Rth and Rr are shown in 
Figure ^ After a rapid increase in the beginning, Rth 
decreases to an asymptotic value slowly. This rela- 
tively low thermalization rate in the later part is due 
to a growth of accumulated mass elements in the col- 
lision area which enlarges the effective collision cross 
section and lessens the relative velocity between two 
streams. In case of run 7, Rth reaches as high as 20 % 
in the beginning, then it decreases to an asymptotic 
value of about 12 %. On the other hand, the cool- 
ing fraction Rr of the same simulation has a gentle 
maximum at t « 1.5 x 10* sec and slowly converges to 
^ 0.4 %. From these ratios, we expect that the overall 
amount of energy that escapes from our calculation 
region as radiation during the whole collision phase is 
^ 0.4 % of the total energy input or ~ 3 % of total 
thermal energy. Note that Rr is just the 'inside the 
box' value. The total luminosity from the whole sys- 
tem including the cylindrically expanding, long PCG 



outside the box will be discussed in 54.2 



The overall evolutionary aspects of Rth and Rr for 
other runs are similar to those for run 7 qualitatively, 



and are different only quantitatively. First, run 1, 
for which no radiation is considered and 7 — 5/3 is 
adopted, has significantly less Rth than run 7. This 
is mostly due to the dependence of the shock con- 
ditions on 7 value: from one-dimensional analysis, 
1^7=4/3/^7=5/3 = 64/49 where u is the post-collision 
internal energy density. However, the actual Rth ra- 
tio between two simulations is quite larger than 64/49, 
probably because the different physical environments 
in the collision region initially induced by different 
7 values continuously make the subsequent coUisional 
evolutions of the two simulations differ. 

Runs 8 and 9 have A's different from that of run 
7. As shown in the previous section, the collision of 
run 9 is more of bouncing while that of run 8 is more 
of merging, which provides more effective thermaliza- 
tion environment. This explains relatively higher Rth 
value of run 8. However, interestingly, the Rr values 
of these runs and run 7 are nearly the same despite 
of considerably different structures in the collision re- 
gion. This is because the less effective collision forms 
a less dense PCG, which quickly radiates thermal- 
ized energy due to a smaller optical depth. The rel- 
atively high Rth and low i?,. values of run 10 can be 
attributed to the same reason. On the other hand, 
both Rth and Rr of run 2 converge to the values very 
close to those of run 7. This implies that the stream 
height is a relatively minor factor in determining Rth 
and Rr- Finally, the amounts of both thermalization 
and radiation decrease as the z offset increases. 

4. DISCUSSION 

4.1. Orbital Energy and Angular Momentum 

The collision changes the kinetic energy and angu- 
lar momentum of the streams. For the calculation of 
the latter, one needs the distance from the collision 
region to the black hole, r^h, and the angle between 
the velocity vector of the stream 1 and the vector to- 
ward the black hole from the collision region, Obh- We 
have adopted r;,^ = 6000 i?© = 4.2 x 10^'' cm from 
Kochanek (1994) and 9bh = 16° from the require- 
ment of the angular momentum conservation. These 
parameters depend on the mass of the black hole and 
the collision impact parameter. 

The density-weighted frequency distribution func- 
tions (DFs) of post-colhsion kinetic energy f{Ek,post) 
and those of z component of the angular momentum 
gjJz.post) for runs 7, 8, and 9 are shown in Figure 
id. Since the orbital evolution of the stream mate- 



rial is determined by Ek and Jz values in the post- 
collision expanding phase, only the cells whose dis- 
tances from the collision center are within a certain 
range (30 < r/lceii < 40) and whose velocity vectors 
are heading away from the collision center were con- 
sidered in the calculation of the DFs. Furthermore, 
since the expanding aspects of PCGs A and B are 
considerably different, we have calculated two sepa- 
rate DFs for the material with an angle between -|-y 
axis and velocity vector projected onto the x-y plane, 
(j), smaller than 225 ° (PCG A) and for the material 
with larger than 225° (PCG B). The f{Ek,post) and 
g{Jz,post) of PCG A (B) are normalized such that the 
pre-collision kinetic energy Ek,pre and z component 
angular momentum Jz,pre values of stream 1 (2) be- 
come unity, respectively. 

Generally, f{Ek) of PCG A are narrower and the 
peak is closer to the unity (less affected by the col- 
lision) than those of PCG B. As discussed in § [3.l| , 
this is because stream 1 has greater M/dA. The DF 
peak deviations from the unity of different runs are 
well explained by the morphology of the collision re- 
gions. While most material of PCG A of run 8, which 
provides the strongest collision geometry, has 20-30% 
less Ek than stream 1, only few percent of stream 1 
of run 9, which has the weakest collision geometry, 
has experienced the Ek reduction of only about 10 %. 
In runs 7 and 8, PCG B has gained Ek during the 
strong collision process. On the other hand, in run 
9, most of PCG B has lost Ek rather than gained. 
This is probably because the bouncing phenomena 
with a small angle at the collision region in run 9 
does not effectively transfer the Ek from stream 1 to 
stream 2 and because the longer shock front between 
two streams has preferentially thcrmalized stream 2 
which has much smaller M/dA than that of stream 
1. 

Unlike f{Ek,post), g{Jz,post) of runs 7, 8, and 9 
shows similar aspects with slightly different peak loca- 
tions. Streams initially have the same (positive) sign 
for Jz, but the collision produces negative Jj, result- 
ing an increase of the total absolute angular momen- 
tum. In general, PCG has 2-3 times greater \Jz\ than 
the initial values Jz,i- Hence the collision changes Jz 
of the stream much more easily than it changes Ek, 
which confirms Kochanek's (1994) estimate. This rel- 
atively large change in Jz is due to the rapid growth 
of the sine function (in mrbhVSuiObh) at small angle 
(Obh — 16°). The average and the standard deviation 
of Jz,post of PCG A in run7 are 1.73 Jza and 0.87 Jza, 



respectively, where only the gas elements whose sign 
of J z, post is the same as Jz,i were considered. Note 
that all the gas elements in our simulations are bound 
to the black hole even after a part of them gains some 
kinetic energy through the collision. This is because 
the simulations here are for the very beginning part of 
the debris stream, which is already well bound before 
the collision due to 1) its proximity to the black hole 
during disruption and 2) hydrodynamic effects that 
convert the orbital energy of the stream into internal 
energy during the revolution around the black hole. 

The exact details of the circularization process is 
still uncertain. If the angular momentum is con- 
served during the circularization process, the process 
will form the debris into a torus with a radius of 
jpost/GMbh, where jpost is the angular momentum 
per mass of the PCG. Then the time scale for that 
process will depend on the difference between the ini- 
tial orbital (kinetic -I- potential) energy per mass of 
the PCG and the energy per mass in a circular orbit 
of the same angular momentum, Ae: 

1 fGAhh^^ 



^post,i 



+ 



^ \ Jpo 



(7) 



where epost,i is the initial orbital energy per mass of 
the PCG. The stream-stream collision lessens Ae by 
changing their pre-collision kinetic energy into ther- 
mal energy (thus decreasing orbital energy) and in- 
creasing angular momentum. We find that for runs 7, 
8, and 9, the collision lessens Ae by a factor of 3 to 10, 
compared to the values which the stream would have 
without the collision. Thus if the angular momentum 
is conserved during circularization process, the time 
scale for that process will be shortened significantly by 
the collision, but the quantitative estimation of such 
decrease in the time scale is still a difficult problem. 

4.2. Luminosity Evolution 

The luminosity from the gas is simply given by 
sum of the cooling rate, equation (||), over the whole 
volume. The luminosity from the gas within the sim- 
ulation box Li,ox{t) for run 7 is plotted in Figure |l2|(a) 
with a solid line. This curve differs from Rr of Fig- 
ure ^ in a sense that the latter is time-cumulative 
energy ratio inside the box. After a steep rise up 
to ~ 10"*^ erg sec in ~ 10^ sec, the luminosity de- 
creases to the steady state value of ^ 10^° erg sec. 
However, since the PCG continues to radiate beyond 
the simulation box, the radiative cooling outside the 
box should also contribute to the total luminosity. 



To describe the behaviour of the gas outside the 
simulation box, we use simple physical arguments. 
The behaviour of the PCG outside the box will be de- 
termined by the competition between the energy loss 
by radiation and that by adiabatic expansion. Sup- 
pose adiabatic expansion is dominant over radiative 
cooling. Since the outgoing PCG outside the box is 
not bounded by incoming streams or shocked gas, it 
may be described by the simple pressure driven ex- 
pansion of a homogeneous gas cloud: 



-1/7 



dP 



-powdwj 



(8) 



where P and p are the pressure and the density inside 
the PCG, w is the time derivative of the width of the 
PCG slice, w, and subscript is for the initial values. 
Due to the one-dimensional nature of the PCG, we 
expect the expansion to be cylindrical: 






(9) 



Combining equations (||) and (J^), one obtaines the 
equation of motion for w: 



dw 
'dt 
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(10) 



Here, the initial width of the PCG slice is defined by 
wq = (ttpowpcg/Mpcg)^^^^, (11) 

where vpcG and Mpcc are the stream velocity and 
the mass flux of the PCG. 

Numerical integration of equation ( |lO| ) with the 
initial values from the results of run 7 is plotted in 
Figure ^ (sohd hue). The width w{t) grows slowly 
in the beginning, but soon attains the constant ex- 
pansion velocity which is the consequence of equa- 
tion (H) . Also plotted are the evolution of the inter- 
nal energy in the slice by pure adiabatic expansion 
u{t)w{t)'^ (dashed line) and by pure radiative coohng 
Uoexp(— t/td) (dotted line), where td oc TesW oc nw^ 
is constant for cylindrical expansion. The former de- 
creases much faster than the latter. Therefore, the 
internal energy will be mainly spent to the adiabatic 
expansion and the PCG may well be described by the 
pressure driven adiabatic expansion. 

It takes approximately tdyn — 5 x 10^ sec for the 
PCG to reach the black hole from the collision re- 
gion. We found from the above w{t) calculation that 



7e// remains very close to 4/3 at least for tdym and 
thus radiation energy will dominate over gas energy 
during that period. When the PCG passes near the 
black hole, it could be significantly distorted by the 
tidal force and adiabatic expansion may not be a good 
approximation because the tidal effect could heat or 
cool the PCG if it passes very close to the hole. How- 
ever, according to the w{t) calculation, the PCG by 
then would have cooled down to T ^ 3 x 10** K. At 
this temperature line opacities will dominate over the 
continuum ones and the PCG will be very optically 
thick and radiates very little. Since we already as- 
sumed that the pre-coUision gas stream maintained 
near this temperature does not radiate, we assume 
a slice of the PCG is luminous only for tdyn after it 
leaves the collision region and ignore its contribution 
to the total luminosity afterwards. 

The total luminosity, L, is the sum of the radiation 
from the gas within the simulation box, Ljyox^ and that 
from the gas outside the box. For the emission from 
the PCG outside the box, we use the results from the 
w{t) calculation. The emission rate per unit length 
slice along the PCG, IsUce, is irw^e. Since TesW is 
constant for cylindrical expansion, the emission rate 
becomes 

hucei-s) (X u{At)w{Atf , (12) 

where s is the distance along the PCG from the 
boundary of the simulation box. At is the time af- 
ter the PCG leaves the boundary, and s = vpca^t- 
Proportional constant is fixed to match the condition 
at the simulation boundary, and the luminosity from 
the PCG outside the box is now the sum of the slice 
emission along the whole PCG. Thus the total lumi- 
nosity of the system is 



m 



-'box 



(t) 



lsliceis)ds. (13) 



This total luminosity is plotted in Figure [l2|(a) 
with a dashed line. While the luminosity inside the 
simulation box (solid line) converges to an asymp- 
totic value, the total luminosity grows as the length 
of the thermalizcd PCG increases. We may extrap- 
olate the total luminosity evolution by assuming a 
steady state for the gas inside the simulation box for 
t > 3 X 10^ sec. Figure [l^(b) shows such extrapola- 
tion, where the total luminosity reaches at t ~ tdyn 
a steady state value of ~ 10"*^ ergsec"^. Since tdyn 
is about a half of the collision phase, the durations 
of the exponential luminosity increase and the steady 
state will be about the same. 



The collision between the incoming stream and 
largely expanded PCG after its passing the black hole 
for the second time will be highly ineffective. This 
ineffective collision phase will last until the stream, 
slightly perturbed by such an inefficient collision, 
passes the black hole and collides with its tail at far- 
ther distance from black hole than the previous colli- 
sion. Since the velocity dispersion of the stream while 
passing the black hole results in a larger cross section 
for more distant collision, subsequent collisions will 
be less and less energetic. The luminosity evolution 
with a shape of Figure |l^(b) will be repeated inter- 
mittently, but with decreasing intensities. After the 
debris is circularized, the viscous accretion process 
will determine the evolution of the luminosity from 
the system. 

Emitted photons from the PCG will look Planckian 
because the gas is effectively optically thick. Since the 
temperature is less sensitive to the expansion of the 
PCG (T oc n^/'^), the characteristic temperature of 
the whole PCG will not be far from that inside the 
simulation box (order of 10^ K). However, it is certain 
that this characteristic temperature decreases as the 
length of the PCG increases. 

5. SUMMARY 

We have performed simulations of supersonic colli- 
sions between two gas streams produced through the 
tidal disruption of a normal star by a super massive 
black hole with a three-dimensional hydrodynamics 
code based on the TVD scheme. A simple treatment 
for radiative cooling valid in any optical depth was 
included in the calculation. 

Thermalized gas emerges out of the narrow, slab- 
like collision region and expands out into a larger 
volume in two opposite directions, forming two ex- 
panding cylinders several times thicker than the pre- 
collision streams. The angle of the shock slab, which 
determines the effectiveness of the collision, depends 
primarily on the mass ffuxes per unit area of the col- 
liding streams. When there is no Lense-Thirring pre- 
cession, approximately 10% of the initial kinetic en- 
ergy of the streams is converted into thermal energy 
during the collision depending on the stream parame- 
ters, and 3 to 4% of the thermal energy is immediately 
radiated away at the collision area in form of radia- 
tion. Off-centered collision which may be caused by 
Lense-Thirring precession can somewhat reduce the 
amount of thermalization and radiation, but still sig- 



nificantly alters the angular momentum vector of the 
outgoing stream. 

The collision alters the angular momentum of the 
stream more easily than it changes the kinetic energy. 
The angular momentum of the incoming stream is 
increased by a factor of 2 to 3, and such increase, 
together with the decrease in kinetic energy, reduces 
Ae, the difference between the initial orbital energy 
per mass of the PCG and the energy per mass in a 
circular orbit of the same angular momentum, by a 
factor of ~ 10. 

The luminosity evolution of the system has three 
phases: the initial burst phase due to effective col- 
lision in the beginning, exponential increase phase 
when the length of the thermalized, adiabatically ex- 
panding post-collision gas increases, and steady state 
phase that lasts until the end of the collision phase. 
The radiation from the system is expected to be close 
to Planckian with effective temperature of ^ 10^ K. 
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y^code) 


7 
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5 


3 


0. 


4/3 


run 9 
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12 
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run 2 
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run 10 


10 
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10 
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run 11 
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4/3 


run 12 


5 
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5 


3 


3.2 
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0. 
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Fig. 1. — Normalized, density- weighted frequency 
distribution of post-collision •y^ff for run 7 at i = 
4 X 10^ sec. Most PCG has jeff very close to 4/3, 
the adiabatic index value for radiation. 
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Fig. 2. — Density, temperature, velocity vector plot 
of z = 50 plane for run 7 at t = 4 x 10'* sec. Units 
of x, y, and z are the grid size of the simulation, 
5.93 X 10^° cm. The contour is for density, grey-scale 
map for temperature, and arrows for velocity vectors 
projected onto the x-y plane. The temperature scale 
is shown as a bar at the top of the Figure. Stream 

1 moves from top to bottom at x = 50 and stream 

2 moves from lower-right corner to upper-left corner. 
Two streams are colliding at the center of the Figure, 
and thermalized PCG is expanding out mainly toward 
lower-left corner. Profiles of the line connecting points 
Pi and P2 are shown in Figure |5[ 
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Fig. 3. — Number density (n), radiative cooling per cell (A), temperature (T), total optical depth (rtot), and effective 
optical depth (r*) profiles of z=50, y=f 5 (a), y=45 (b), y=55 (c), and y=60 (d) lines for run 7 at t = 4 x 10^ sec. 
The units for n, A, and T arc cm^'^, fO^^erg/sec, and K, respectively. 
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Fig. 4. — Cooling map in grey-scale for z = 50 plane 
of run 7 at t = 4 X 10'* sec. Overplotted is the density 
contour, and a grey-scale bar is shown at the top in 
units of erg /sec. 
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Fig. 5. — Number density (n), total internal energy 
per cell divided by a time step {Eth/ /S.t), cooling per 
cell (A), temperature (T), total optical depth {Ttot)i 
and effective optical depth (r*) profiles along the line 
connecting points Pi and P2 in Figure 0. The units 
are same as in Figure with Eth/^t having the same 
units as A. 
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Fig. 6. — Density, temperature, velocity vector plot 
of z = 50 plane for run 8 at i = 4 x 10* sec. See 
Figure g for explanations. 



Fig. 7. — Density, temperature, velocity vector plot 
of z = 50 plane for run 9 at t = 4 x 10* sec. See 
Figure y for explanations. 
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Fig. 8. — Three-dimensional isodensity surface plot 
for run 12 at i = 4 x 10'' sec. The surface represents 
n ~ 8 X 10^'''cm^'' points. The view angle is from +x 
axis to —X axis. Stream 1 moves from right to left in 
z = 50 plane, and stream 2 moves from left to right 
with an angle of 45° into the Figure in z = 58 plane. 
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Fig. 9. — Evolution of the ratio of total thermal en- 
ergy to total input kinetic energy {Rth', solid lines) 
and the ratio of total cooling inside the simulation 
box to total input kinetic energy {Rr', dashed lines) 
of all simulations. 
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Fig. 10. — Normalized, density- weighted frequency distribution of post-collision kinetic energy Ek,post (a), and 
z-component post-collision angular momentum Jz,post (b). Thick lines are for the gas with an angle between -|-y 
axis and its velocity vector projected onto x-y plane, 0, smaller than 225° (PCG A), thin lines for the gas with (p 
larger than 225° (PCG B). Ek.post and Jz,post of PCG A (B) are in units of pre-collision kinetic energy of stream 1 



(2), Ek.pre, and pre-collision angular momentum of stream 1 (2), 
of PCGs A and B are normalized separately. 



Jz,pre, respectively. The frequency distributions 
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Fig. 11. — Evolution of the width of the PCG shce ■w{t) (sohd). The PCG shcc is assumed to be a homogenous gas 
cloud that expands adiabatically and cylindrically. Initial conditions are obtained from the results of run 7. Also 
plotted are the internal energy evolution of the PCG slice by pure adiabatic expansion, oc u{t)w{t)'^ (dashed), and 
that by pure radiative cooling, oc exp(— t/t^) (dotted). All variables are plotted in logarithmic scale and normlized 
to their initial values. 
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Fig. 12. — (a) Temporal evolution of the luminosity inside the simulation box (solid line) and the total system 
luminosity including the estimation for the PCG (dashes line) for run 7. (b) Extrapolation of total luminosity 
evolution with the physical quantities being the steady state values at the last phase of the simulation. A slice of 
the PCG is assumed to adiabatically expand in two dimensions and to be luminous only for tdyn after it leaves the 
collision region. 
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